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ABSTRACT 

Gamma-ray catalogs contain a considerable amount of unidentified sources. Many of these are 
located out of the Galactic plane and therefore may have extragalactic origin. Here we assume 
that the formation of massive black holes in galactic nuclei proceeds through a quasi-star stage and 
consider the possibility of jet production by such objects. Those jets would be the sources of collimated 
synchrotron and Compton emission, extending from radio to gamma rays. The expected lifetimes of 
quasi-stars are of the order of million of years while the jet luminosities, somewhat smaller than that 
of quasar jets, are sufficient to account for the unidentified gamma-ray sources. The jet emission 
dominates over the thermal emission of a quasi-star in all energy bands, except when the jet is not 
directed towards an observer. The predicted synchrotron emission peaks in the IR band, with the 
flux close to the limits of the available IR all sky surveys. The ratio of the 7-ray flux to the IR flux 
is found to be very large (~ 60), much larger than in BL Lac objects but reached by some radio- loud 
quasars. On the other hand, radio- loud quasars show broad emission lines while no such lines are 
expected from quasi-stars. Therefore the differentiation between various scenarios accounting for the 
unidentified gamma-ray sources will be possible at the basis of the photometry and spectroscopy of 
the IR/optical counterparts. 

Subject headings: black hole physics; accretion; galaxies: active; gamma ray bursts 



1. INTRODUCTION 

A significant number of sources in the 7-ray 

sky has no counterp arts at other energy bands 

(iCasandiian fc Grenier I [20081: lAbdo et al. 1 120101: 



0! 



Bir d et al. II2010T) des pite considerable efforts to identify 
them (|Maeda et al. 1120111 iMasetti et al. 112010ft . This 
puzzle has persisted in gamma-ray astronomy since 
many years. Even the source catalog of the Fermi-LAT 
instrument with the best position errors contains 
about 300 of stil l unclassified, high-latitude sources 
(|Abdo et al. 1 120101 . The nature of these sources is thus 
difficult to establish but because of their positions they 
are most likely of extragalactic origin. This means that 
explaining their nature requires the identification either 
with known active galaxies or the existence of a large 
population of a new type of extragalactic sources. The 
first option has been recently considered by Massaro 
et al. (2012a; 2012b). We address here the second 
possibility. 

The part of the galaxy evolution which still remains the 
most mysterious is the pre-quasar epoch of formation of 
a massive black hole. The process must be very rapid 
since massive black holes at galactic centers must form 
surprisingly early in order to ac count for the quasar-type 
activity at ver y large redshifts ()Volonteri fc Reesl 120061 : 
IHaimanl [2010^1 . More and more distant quasars located 
at very high redshifts are continuously being discov- 
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ered, and their central black hole masses are large (e.g 
ULAS J112001.48+064 124.3: z = 7.085, M = 2 x 10 9 M, 
(|Mortlock et al. Il20ll j. 

An attractive scenario for the growth of black holes 
in time considerably sh orter than the Sa l peter timescale 
has been proposed by [Begelman et al. I (|2006f ) . There 
the early accretion proceeds through a quasi-star stage, 
with a massive accreting envelope surrounding an ini- 
tially small, few solar mass seed black hole. The 
detailed structure and quasi-stationary evolution of a 
quasi-star, assuming spherical sym metry, has been re - 
cently discussed in more det ai l by IBegelman I ( 2010ft , 
I Volonteri fc Begelman I (pOTOft . IBall et al. I (|2011l ). and 
iBall et al. I (j2012ft . However, as n the case of a collaps- 
ing star, the configuration does not have to be spher- 
ically symmetric. If the quasi-star material posseses 
some angular momentum, a jet may form, as in gamma- 
ray bursts or active galactic nuclei. A strong collima- 
tion of radiation may lead to enhancement of the ob- 
served radiation, large enough to make 7-rays detectable 
from cosmological distances and to make domination of 
IR/optical synchrotron radiation over uncollimated ther- 
mal radiation from the quasi-star's envelope. Finding 
such IR/optical counterparts to the 7-ray sources will 
provide oportunity to justify the model. In the present 
paper we consider the possibility of the quasi-star jet for- 
mation and its observational consequences. 



2. MODEL OF A QUASI STAR 
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The formation of seeds for the present-day massive 
black holes is unclear and a number of possi bilities 
are un der discussion since the seminal paper by iRees I 
(1978), and various mechanisms may actually coexist. 
At la rge redshifts (z > 20) P opulation III stars likely 
form (jBromm fc Larson"ll2004h and finish their evolution 
as black hole s with masses of o rder of a few hundred 
solar masses (jHeger et al. 1 120031 ). At intermediate red- 
shifts, where young galaxies form, protogalaxy mergers 
lead to creation of a strong gas inflow toward the cen- 
ter and a massive quasi-star may form. Such a spher- 
ical configuration is unstable but if partially supported 
by rotation, may form a quasi-stationary configuration 
of a quasi-star. A black hole at the center is accreting 
at highly super-Eddington rate and providing the en- 
ergy support for the expanding self-gravitating envelope; 
the structure is therefore similar to that of a red giant 
star but with much higher mass. The inner black hole 
slowly grows till it contains about 1-10 per cent of the 
initial quasi-star's mass when the configuration looses 
stability and becomes dispersed, leaving behind an in- 
termediat e mass black hol e . Suc h a scenario was pro- 
posed bv iBegelman et al~l (l2006|) and subsequently de- 
veloped in a number of pap ers iBegelman et al. 1 ,2008; 
iVolonteri fc Begelman 1 12010: 1 Ball et al. II2011D . " 

In the present paper we supplement the quasi-star 
scenario with the possibility of a jet formation which 
will shape its observational appearance. In this case 
the jet production is presumably similar to that occur- 
ring in AGN and GRBs, i.e. mediated by rotation of 
poloidal magnetic fields in the BH and/or disk magne- 
tosphere. Such fields can effectively be transported from 
the outer regions to the center only by geometrically 
thick accretion flows flCao 11201 lHMcKinnev et al. 1120121 : 
iTchekhovskov et al. Il2012t ). 

Here we use a very simple parametric description of a 
quasi-star and the jet and study the observational con- 
sequences of this model. 

2.1. The quasi-star structure 

A quasi-star (hereafter QS) in our scenario consists of 
two parts: a spherically symmetric envelope radiating 
thermally at the Eddington luminosity and a relativistic 
jet. Such a structure can be described using very few 
parameters. The schematic view of our model is shown 
in Fig. m 

2.2. The spherical component 

The basic parameter of the object is the quasi-star 
mass, A/qs- We adopt as representative the value of 
10 7 Af Q and further use the total mass in the dimension- 
less units ofMQgj. 

The QS isotropic thermal luminosity is simply deter- 
mined by the condition of the Eddington luminosity at 
the surface of the spherical envelope, 

L QS = 1.3 x 10 45 M QSi7 [ergs" 1 ]. (1) 

The value of the central black hole mass, the surface tem- 
perature and the radius of the spherical configuration 
evolve. The mass of the central black hole grows, the en- 
velope expands and the effective temperature drops. In 
general, the temperature is determined by the QS mass 
and the central bla ck hole mass through a simp le for- 
mula (see Eq. 4 in IVolonteri fc Begelman I I2010D which 




QS 



Fig. 1. — The structure of a quasistar. Representative parame- 
ters: mass ~ 10 7 Mq, luminosity ~ 10 45 erg s _1 , black hole mass 
10 5 Mq, and the geometrically thick, advective disk with accretion 
rate ~ 2A/Qyr — 1 , with radiative efficiency of 0.01. Accretion en- 
ergy is converted to power the bipolar jets with an efficiency of 
0.1 and the jets have Lorentz factor T ~ 15, radiative efficiency 

0. 1 and an isotropic luminosity ~ 10 48 erg s — 1 . The jets produce 
gamma rays in the reconfinement shocks formed within 0.01-1 tqs, 

1. e. 10 15 - 10 17 cm. 



results from the polytropic approximation to the stellar 
structure, and the configuration is described by two pa- 
rameters. In order to simplify the parameterization we 
use the final stage of the QS as representative. This final 
stage is achieved when the star reaches the Hayashi limit, 
roughly given by the effective temperature TQs, m in, equal 
to 3000 K for a solar metallicity. We can then determine 
the black hole mass of s uch c onfiguration from Eq. 4 in 
IVolonteri fc Begelman I tf2010h . 



BH 



.5 X 



iU M QS,7 e d,-2 



(2) 



where £d is the radiative efficiency of the disk. We ex- 
pect this efficiency to be low because the accretion rate 
onto the central black hole is highly super-Eddington and 
the flow is advection-dominated, so we take the value of 
0.01 as representative. The corresponding accretion rate 
(£qs = e d M d c 2 ) is 



Lqs 
e d c 2 



2.2M QS!7 £ d 1 _ 2 [Afoyr- 1 ] 



(3) 



i.e. the Eddington ratio is of order of 100 at the end of 
the QS evolution, and much higher at earlier stages. 

The thermal emission of the QS should be visible in 
the IR, in 2 to 3 micron wavelength range (K and L 
bands), taking into account the red color of the envelope, 
and the redshift. The assumption of TQs.min = 3000if 
determines also the radius of the QS 



TQS 



Lqs 



47rCTT QS,min 



1.5 x 10 17 M^ 7 [cm] 



(4) 



where a is Stephan-Boltzman constant. The object is 
very extended due to the rapid decline in the density 
with radius. 
The configuration is long-living: 
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although it lasts less than the quasar activity period. 
2.3. The jet 

We further assume that a fraction of energy of the ac- 
creting material is converted to power the relativistic jet 
with an efficiency of r]j . For a super-Eddington accretion 
flow the radiative efficiency e<± is diminished by advection 
but there are no such constraints for a jet, and rjj can be 
higher than so we take the value of 0.1 as a scaling 
factor. Therefore the total power of a quasi-star jet can 
be parameterized as 



P j =r ?j M d c 2 = ^L QS 
= 1.3 x 10 46 Af QS ,7 



- [erg s 2 ] 



(6) 



For such a jet power the breakout condition, tQs ^> 
r Qs/ v h, is very well satisfied, where Vh is the velocity 
of the jet head paving the way outwards through the 
QS and determined by the balance between momenta 
fluxes of a jet and external medium as measured in the 
head comoving frame (see, e.g. Scheuer 1974; Waxman 
& Meszaros 2003). The head velocity is relatively slow 
when emerging from the star, Vh ~ 430 km s _1 , but it is 
still highly supersonic due to the low temperature of the 
envelope. 

If the radiative efficiency of the jet is ej and its radi- 
ation is collimated into a solid angle f2 ra d = 7r0^ ad , then 
an observer located within this angle will detect a flux 
F = L^iso/And^, where the isotropic jet luminosity is 



4tt 



1 £ r . 



ad 

1.1 x 10 48 - 1 



Ed, 



(0rad/4°) 2 



[erg s 



(7) 



Here 4° is the maximum angle at which the condition of 
the causal connection between the two opposite sides of 
a jet is satisfied for a jet Lorentz factor r ~ 15 - a typical 
value for jets in blazars associated with quasars. Hence, 
the observed flux from the jet is by a factor 



L 



1.1SO or 1 a2 —1 

- 8.5 x 10"e j ,_i7^_ie d i _ 2 



QS 



(flrad^T 
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larger than the flux from the QS. The observed jet emis- 
sion always dominates over the thermal emission of the 
spherical envelope if the jet is directed towards an ob- 
server. 

Accretion disks formed in the center of QSs are 
super-Eddington (see Eq. |3] and the discussion below) 
and therefore geometrically thick. Such disks can ef- 
ficiently advect poloidal magn etic fields t o the central 
regions of the accretion flow (|Cao 1 120111 ). This may 
lead to formation of powerful, Poynting flux dominate d 
jets (jMcKinnev et al. ll2012tlTchekhovskov et al. 1120121 ). 
Such jets undergo acceleration up to the moment when 
about half of their electro-magnetic power is converted 
to the kinetic energy of the carried matter, i.e. when 
cb = Lb/Lk drops to unity, where Lb is the mag- 
netic energy flux, Lk is the kinetic energy flux, and 
Pj = Lb + Lk- This happens at a distance r ~ 10 2 — 
10 3 r g , where r g = GM B h/c 2 (jKomissarov et al. I [20071: 



iTchekhovskov et al. I [200l iLvubarskv I [2010h . Further 
away the conversion proceeds much slower, but possibly 
at r ~ 10 4 — 10 5 r g the ratio <jb becomes small enough 
to allow for the form ation of strong shoc ks and efficient 
particle acceleration (jSikora et al. 1 120051 ). Because the 
size of our QS at t = £qs is given by Eq. [4] while the 
gravitational radius of a black hole with the mass given 
by Eq. 2 is 



1.3 x 10 



li 



M 7/8 -1/2 



cm] 



(9) 



we have a very large ratio of the envelope to the black 
hole radii rQs/V g ~ 10 6 . Therefore one may expect that 
efficient shock formation and non-thermal radiation pro- 
duction takes place within a distance range 0.01 — 1 tqs, 
where due to the interaction with the QS's envelope the 
reconfinement shocks develop. 

3. RADIATIVE PROPERTIES 

We now estimate the properties of radiation emerging 
from a QS with a jet in order to formulate the observa- 
tional characteristics of such objects. 

3.1. Radiative processes 

The dominant radiative processes are: synchrotron, 
synchrotron-self-Compton (SSC), and external- 
radiation-Compton (ERC). Radiative electron energy 
losses in these processes scale with the magnetic energy 
density, u' B , synchrotron energy density, u' syn , and the 

external radiation energy density, u' ext , respectively, 
all measured in the jet co-moving frame. Noting that 
cu' B , Kr 2 {9- i T) 2 and that for <tb <C 1, Lb — fB-Pj, 



we have 



and 



r MB** = -■ 



4 rL Qs r 



(10) 



(11) 
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we find the ratio of the external to the synchrotron flux 
equal to 



Perc £erc e d ,-2T(r/15) 2 (0jr) 2 
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syn 



CB,-l»7j,-l 



(12) 



where #j is the opening angle of a jet. The parameter r 
scales the real density of QS radiation entering a jet at a 
distance r by a factor LQg/47rr 2 c. At r = tqs it is r ~ 1, 
and at r < tqs it is t > 1. Its dependence on a distance 
within the QS envelope is determined by the dependence 
of the opacity on radius as well as on multi-scatterings 
of QS photons off the jet walls. Making the comparison 
of u' B with u' syn one can find that L syn 3> issc- 

3.2. Radiation spectra 

In the presence of protons, the electrons are ener- 
gized in at least two steps, the preheating and the 
stochastic. In the first step, due to coupling via col- 
lective processes in magnetized plasma with protons, 
electrons are pr eheated up to thermal energies of the 
shock ed protons (|Spitkovskv 11200 8: Sironi & Spit kovskvl 
120111 ). In the second step, electrons together with 
proto ns participate in the Fermi 1 acceleration pro- 
cess (|Blandford fc Eichlerl 119871 : iNiemiec & Ostrowskil 
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2004). For relativistic shocks (Lorentz factor of the up- 
stream flow relative to the shock front is > 2), the elec- 
trons reach ultra-relativistic energies already in the first 
step. This energy in m e c 2 units is 7 m = rem p /m e , where 
for mildly relativistic shocks K ~ 1. Electrons with this 
energy are responsible for the luminosity peak in the vL v 
representation. 

Electrons with such energies produce ERC spectra 
peaking at photon energies 

h^ERC ~ (r 7m ) 2 hi/erf (1 + z)- 1 , (13) 

where h is Planck constant and z is the source cosmolog- 
ical redshift. 

Noting that the peak of the thermal spectrum in the 
vL v representation is f cx t — 3.9 kT, and that 

T = 3000r 1 / 4 (r QS /r) 1 / 2 [K] , (14) 

one finds 

h^ERC ~ 100(r/15) 2 (r QS /r) 1 / 2 r 1 / 4 [MeV] . (15) 

Hence for any location of radiative processes within a dis- 
tance range 0.01 — 1 ?"qs, the dominant radiative output 
is predicted to be within the Fermi energy band. 

The 7 m electrons will produce synchrotron spectra at 
i/ syn ~ 3.7 x 10 6 7 m ST/(l + z) Hz, and for u' B given by 
Eq.[l0] 

B> = 3.94f_ t {r Q s/r) ^ > ( 16 ) 

the peak of the synchrotron radiation produced at z = 3 
is located at frequencies 

Vsyn ~ 2.0 x 10 14 ^— 1/2 t Hz ] ■ 

(17) 

Our results show that QS jets produce synchrotron 
radiation in the IR band and that this radiation dom- 
inates over the QS thermal spectrum. Hence, the di- 
agnosis based on I R sp ectral properties developed by 
I D'Abrusco et al~l (|2012| ) might not distinguish between 
blazars and QSs. However, one may discriminate be- 
tween these two classes of objects noting that we pre- 
dict a ratio Lp CTm i/ Lm ~ several tens. This is observed 
in some luminous blazars associated with quasars, but 
practically never in BL Lac objects. On the other hand 
we can discriminate QSs from quasar-associated blazars 
by noting that in the latter the broad emission lines are 
seen, whereas they are not expected in QSs. 

3.3. Statistics 

The collimation of the jet radiation increases the ob- 
served luminosity but the high collimation of the quasi- 
star jet lowers the probability of its detection. On the 
other hand, all galactic centers should go through this 
phase while only a small fraction of black holes pass 
through a stage of a bright quasar. Therefore statis- 
tically the number of sources is significant and should 
match the fraction of the still unidentified gamma ray 
sources. Below we consider this issue quantitatively. 

The energy flux for the weak unidentified Fermi sources 
is of the order of 10~ n erg cm~ 2 s _1 and such fluxes cor- 
respond to absolute luminosities of 1.1 x 10 48 erg s _1 (see 



eq. [7J at the redshift z = 3.4, in the standard cosmology 
(h = 71, n m = 0.27, n A = 0.73). 

The number of quasi-stars forming at redshifts be- 
tween 2 and 4 has been estimated to be of the or- 
der of 0.01 per cubic Mpc i n the comoving frame by 
iVolonteri fc Begelman I (|2010t ) (see their Fig. 5). With 
this rate we expect a total number of 10 10 quasi-stars 
to form in this redshift range. The number of sources 
to be observed at a quasi-star stage will be much smaller 
since this stage lasts thousands to millions of years which 
is short compared to the age of the universe. A rough 
estimate of the QS phase duration is given by tQg (see 
Eq.lSJ), thus reducing the number of active quasi-stars to 
a representative value of 3.1 x 10 5 for all dimensionless 
parameters equal 1. Since the emission is collimated, the 
number of sources with jets toward an observer will be 
further reduced. In general, combining all coefficients we 
obtain 

N = m(9/^fM-^ J S e 1 d %, (18) 

which shows that both the typical 7 -ray fluxes and the 
number of the high-latitude, uni dentified Fermi sou rces 
in the First Fermi LAT Catalog ()Abdo et al. 112010ft can 
be reproduced in the QS-jet scenario. 

Since the number of on-axis QSs at still larger redshifts 
(from 4 to 10) is only by a factor of 2 larger than the num- 
ber of resolved sources, their contribution to the stil l un- 
explained 7 -ray background (|Inoue fc Totani l[2009h can 
be significant (up to 50%) taking into consideration their 
bolometric luminosity, but more careful analysis should 
be done taking into account the source and the back- 
ground spectral shape. Off-axis sources can contribute 
to the hard X-ray background, but not significantly (a 
few per cent or less). 

4. CONCLUSIONS 

A large fraction of unidentified 7-ray sources is located 
at high latitudes and therefore most of them are consid- 
ered to be extragalactic. We propose that quasi-star jets 
can account for these sources. In analogy to blazars they 
are predicted to produce nonthermal spectra, at lower en- 
ergies (optical-IR) dominated by the synchrotron mecha- 
nism and in the 7 -ray band by the inverse-Compton pro- 
cess. However, they can be distinguished from blazars 
by the following properties: 

(i) the expected ratio of the gamma-ray to the IR com- 
ponents is close to 60, much larger than in BL Lac objects 
(see e.g. Fig. 8 of D'Abrusco et al. 2012) but achievable 
by some luminous blazars associated with quasars; 

(ii) no broad emission lines are expected as in most of 
BL Lac objects but unlike in radio-loud quasars. 

Since the typical Fermi flux of the unidentified sources 
is of the order of 1Q — 11 erg cm~ 2 s _1 , the predicted 
IR/optical magnitudes of the QS counterpart are ~ 
17 — 18 mag in K band and 20 — 21 mag in R band. 
The currently available IR all sky survey WISE has the 
limiting magnitude of 16 — 17, which is at the border 
of the QS detection. Optical SDSS survey also covers a 
significant fraction of th e sky and contains so urces down 
to 22.5 mag in R band (jAihara et al. 1 120111 ). Since the 
IR/optical emission is non-thermal and likely highly vari- 
able, color diagrams and variability can help to select 
proper candidates. 
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